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The  methodology  of  multithreshold  analysis,  a new  approach  to  laser-damage 
research,  is  described.  The  comparison  of  thresholds  for  va-  >ous  damage-related 
effects  identifies  dominant  failure  mechanisms  and  ^ro'  ides  better  guidance  for 
laser-materials  technology.;  After  a brief  description  of  apparatus,  the  proce- 
dure for  routinely  measuring  up  to  eight  different  thresholds  per  sample  is 
given.  The  maximum-likelihood  principle  is  used  to  derive  an  algorithm  for 
computing  thresholds  and  standard  deviations.  The  use  of  a standard  gold  sample 
to  verify  reproducibility  and  to  maintain  long-term  calibration  is  discussed. 
Examples  of  multithreshold  results  on  uncodted  and  coated  infrared  optical 
components  are  presented.-  The  following  are  some  of  the  effects  for  which 
thresholds  are  compared:  slip,  roughening,  cracking,  pics,  melting,  craters, 
delamination  of  coatings,  ion  and  light  emission,  and  work  function  change. 

Key  words:  Damage  threshold'.;  failure  modes;  laser-optical  components;  Mo 
mirrors;  optical  coatings;  pitting;  ripples;  standard  sample  Au  mirrors;  sur- 
face characterization;  surface  uniformity;  threshold  estimation;  windows. 


1.  Introduction 

The  damage  resistance  of  optics  for  lasers  may  be  Improved  with  an  understanding  of  Important 
damage  mechanisms,  and  their  relationship  to  material  characteristics.  The  comparison  of  thresholds  for 
the  various  observable  damage-related  effects  provides  a good  basis  for  this  understanding..  Some  of 
the  effects  typical  of  a bare  metal  mirror  may  be  seen  in  figure  1.  Slip  bands  and  intergranular  slip 
indicate  susceptibility  to  stress  damage,  which  can  be  related  to  yield  strength,  as  shown  in  a compan- 
ion paper  [1]^.  Melting,  which  produces  the  dark  area  in  the  dark-field  micrograph,  is  an  indicator  of 
optical  absorption.  Laser-induced  pitting,  which  occurs  mainly  on  grain  boundaries  in  figure  1,  can  be 
caused  by  segregated  impurities.  Cratering,  identified  by  a raised  rim  around  the  melt  zone,  is  usually 
caused  by  pressure  associated  with  vaporization.  Other  effects  measured  by  other  techniques  may  also  be 
informative.-  A change  in  surface  work  function,  for  example,  can  signal  changes  in  surface  composition 
or  topography..  Ion  and  lipht  emission  indicates  plasma  formation.. 

The  present  work  presents  the  methodology  of  obtaining  thresholds  routinely  and  accurately  for  as 
many  as  eight  different  laser-induced  surface  effects..  After  a brief  description  of  the  apparatus, 
details  of  which  ar'  given  elsewhere  [2],  the  procedure  for  acquiring  and  reducing  the  data  is  discussed.; 
The  derivation  of  an  algorithm  based  on  the  maximum-likelihood  principle,  which  is  used  to  compute 
thresholds  and  standard  deviations,  is  also  given.  We  then  discuss  the  use  of  a standard  gold  sample  to 
verify  reproducibility  and  to  maintain  calibration..  Finally,  we  present  typical  results  and  their 
interpretation  on  various  infrared  optical  components. 

2.  Apparatus 

Measurements  are  made  in  a previously  described  laser  damage  facility  [2],  whe~e  comprehensive 
target  monitoring  and  careful  beam  characterization  are  key  features.  Our  ultrahi„h  vacuum  (UHV)  test 
chamber  was  equipped  with  the  following  monitoring  apparatus  for  the  results  presented  here:  (1)  20- 
power  microscope;  (2)  Faraday  cup  with  grounded  entrance  grid  (10-mA/sr  detection  sensitivity  at  35°  off 
the  laser  beam  axis)  for  sensing  of  target-emitted  ions  and  neutrals  having  sufficient  energy  to  produce 
secondary  electrons;  (3)  an  Auger  analyzer  with  electron  imager  [2],  The  electron  images  are  generated 
by  work  function  variations  ^-0.1  eV  over  the  target  surface  caused  by  variations  in  surface  composition 
and  topography. 

Figure  2 illustrates  the  Gaussian  spatial  intensity  distribution  of  the  10.6  um  laser  beam  when 
focussed  on  the  target  by  a ZnSe  lens  having  a focal  length  of  24  cm.-  Temporally,  each  pulse  consists 
of  a train  of  mode-locked  spikes  which  form  an  envelope  of  100  nsec  nominal  duration.;  Inresholds  given 
are  peak  thresholds,  i.e.,  they  refer  to  the  time-integrated  fluence  of  the  damaging  pulse  at  the  point 
of  maximum  spatial  intensity.; 


* Work  supported  by  the  Naval  Air  Systems  Command,  the  Office  of  Naval  Research  and  NWC  Independent 
Research  Funds. 

1.  Figures  in  brackets  indicate  the  literature  references  at  the  end  of  this  paper.- 


The  data  reduction  equipment  consists  of  a Hewlett-Packard  Model  9821  calculator  interfaced  with  a 
Tektronix  Model  4662  interactive  digital  plotter.  Data  reduction  is  performed  entirely  in  the  calcu- 
lator, while  the  plotter  is  used  only  for  display  of  final  results. 

3.  Experimental  Procedure 


The  sample  is  exposed  to  one  pulse  at  each  of  approximately  30  sites  spaced  at  1-mm  intervals 
along  its  surface.:  Pulse  energies,  which  are  measured  for  each  pulse  [2],  are  uniformly  distributed 
over  a range  which  includes  the  thresholds  for  each  of  the  damage  effects  of  interest.  The  occurrence 
of  light  emission  (flash)  and  pit  formation  are  noted  by  an  observer  as  they  occur,  while  damage  fea- 
tures which  are  more  difficult  to  identify,  e.g.,  slip,  melting  and  cratering,  are  categorized  with  the 
aid  of  more  sophisticated  optical  microscopy  after  the  sample  has  been  removed  from  the  test  chamber 
(fig.  1).,  The  damage  observed  at  eacn  site  is  summarized  by  a set  of  code  numbers,  each  number  repre- 
senting a specific  effect.  The  data  is  reduced  by  entering  the  damage  summary  and  corresponding  peak 
energy  fluence  for  each  site  into  the  calculator  and  applying  the  algorithm  derived  below.. 


3.1  Threshold  Algorithm 

The  routine  determination  of  up  to  eight  thresholds  per  sample  tested  requires  an  algorithm  permit- 
ting rapid  data  reduction.  On  the  other  hand,  effective  use  of  available  data  is  important,  since  the 
fluence  range  of  thresholds  is  often  wide,  and  therefore  incompatible  with  a high  density  of  data  p-  j.nts.; 
In  this  section  we  develop  an  algorithm  which  yields  reproducible  results  with  a minimum  of  computational 
complexity. 


Consider  the  probability  Pj.  of  obtaining  a certain  type  of  damage  times  out  of  shots  on  fresh 
sites  (1-on-l  damage)  [3,4],  all  at  the  same  fluence  This  is  represented  by  the  binomial  distribu- 

tion 


p (X  ;N  ;p  ) » ) P Xi 
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where 


pi  “ 


(2) 


the  probability  of  single-shot  damage  is  a function  of  the  damage  threshold  ♦ and  its  standard  deviation 
o,  in  addition  to  the  fluence  d^.-  The  problem  is  to  obtain  statistical  estimators  for  t and  o in  terms 
of  the  damage  data  Xj_,  Nj_  and  Our  solution  is  based  on  the  method  of  maximum-likelihood  [5).  Maxi- 

mizing the  logarithm  of  the  likelihood  with  respect  to  the  parameters  £ » ( or  t,  leads  to  two  expres- 
sions of  the  fcrm 
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Solving  simultaneously  for  $ and  0 yields  the  corresponding  maximum  likelihood  estimators  i and  a.-  How- 
ever, the  functional  form  of  p^  is  required.-  For  computational  simplicity  we  assume  that  p^  is  uniformly 
distributed,  i.e.  , 


where 


(4) 


(53 


Substituting  eq.  (4)  in  eq.-  (3),  one  finds  that  only  those  terms  J for  which  iu,,  < 1/2  are  non- 
itio.;  To  further  simplify  the  computation  we  neglect  terras  in  the  denominators,  which  are  of  the 
form  i/4  - Uj-.-  As  a result,-  our  estimators  differ  trom  the  true  maximum-likelihood  estimators  in  that 
data  points  near  u^  * 1/2  are  le,-,  heavliv  weighted  vela:ive  to  those  near  uj  - 0.  This  compensates  .to 
some  degree  foj  the  fact  that  the  distribution  given  bv  eq . (■»!  generally  represents  the  sample  distri- 
bution less  accurately  near  m * 1/2  Introducing  e ; . 5)  and  solving  tor  1 yields 


(6) 
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Solving  for  a gives 
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■ 0.,  Noting  that  eq.  (6)  defines  the  sample  mean 


- |«j)  - 0 


(7b) 


It  follows  from  the  definition  of  j that  eqs.  (6)  and  (7)  apply  only  to  data  within  the  interval 


(8) 


While  the  estimators  $ ar  1 a are  expressed  entirely  in  terms  of  the  damage  data,  the  data  interval  given 
by  eq.  (8)  involves  the  unknowns  $ and  a,  and  must  be  determined  by  trial  and  error.-  A satisfactory 
method  for  this  has  been  developec , and  is  Included  in  the  routine  data  reduction  procedure  described 
below. 


3.2  Data  Reduction 

Equations  (6)  through  (8)  are  applied  to  each  damage  effect  in  turn  via  the  calculator  as  follows. 
Flue'.cea  associated  with  the  code  number  corresponding  to  a particular  effect  are  sorted  into  a damag- 
ing, and  the  remaining  fluences  into  a nondamaging  category.  Damaging  fluences  are  then  arranged  and 
listed  in  order  of  ascending  value,  while  the  nondamaging  fluences  are  listed  in  descending  order.  Data 
near  the  head  of  each  list  is  examined  for  consistency  by  the  operator  and  reclassified  or  discarded  as 
required.  Equations  (6)  and  (7)  are  then  applied  to  the  revised  lists.  The  j summations  are  initially 
taken  over  all  data  in  the  smallest  interval  which  (a)  includes  the  lowest  damaging  and  highest  nondamag- 
ing fluence,  and  (b)  Includes  a damaging  and  nondamaging  fluence  as  upper  pd  lower  limits,  respectively. 
The  interval  is  then  extended  as  required  to  provide  the  minimum  value  of  a,  determined  by  applying  eq. 
(7)  to  each  new  interval.:  This  minimum  a and  the  ♦ obtained  by  applying  eq.  (6)  to  the  same 
interval  are  taken  as  the  final  results  for  the  effect  in  question.  In  practice,  it  is  found  that  i is 
rather'  insensitive  to  the  data  interval.  However,  o depends  more  strongly  on  the  number  of  data  points 
used  and  their  separation.  Given  adequate  data  the  major  contribution  to  a usually  comes  from  the  varia- 
tion of  the  damage  threshold  across  the  sample  surface,  with  a relatively  small  contribution  coming  from 
the  experimental  uncertainty.; 

Repeated  application  of  the  algorithm  gives  an  average  threshold,  together  with  its  standard  devia- 
tion, for  each  effect.  The  resulting  thresholds  for  the  target  are  presented  in  the  form  of  a bar  graph 
or  damage  profile,  examples  of  which  are  presented  below. 

3.3  Reproducibility  and  Long-Term  Calibration 

The  reproducibility  of  the  threshold  results  can  be  demonstrated  by  repeated  application  to  a 
sample  with  uniform  damage  characteristics  over  its  surface,  as  indicated  by  relatively  small  standard 
deviations  of  the  damagf.  thresholds.  A mirror  produced  by  electrodeposition  of  Au  on  diamond-turned 
bulk  Cu  has  this  desirable  property.;  The  oxidation  resistance  is  an  added  advantage,  which  makes  this 
sample  especially  attractive  as  a standard  for  calibration  maintenance.  The  ability  to  directly 
standardize  damage  thresholds  is  an  invaluable  asset  in  the  long-term  development  of  damage-resistant 
optical  components. 

Figure  3 presents  the  results  of  two  independent  multithreshold  measurements  on  the  standard  gold 
sample.  Figure  3(a)  was  obtained  before,  and  figure  3(b)  after  replacing  the  output  coupler  on  the 
laser  over  one  month  later.  After  correction  for  a 10X  change  in  focal  spot  diameter,  the  thresholds 
for  all  effects  listed  except  one  agree  within  the  combined  standard  deviations  given.  Statistically, 
approximately  half  the  thresholds  3hould  agree  this  well  when  the  effects  are  independent.-  It  is  clear 
that  changing  the  output  coupler  had  negligible  effect  on  the  calibration..  Furthermore,  we  may  conclude 
that  the  threshold  algorithm  gives  consistent  and  reproducible  results,  the  accuracy  of  which  is 
conservatively  represented  by  tne  standard  deviations. 


3.4  Exaaples 


Muitithresnoid  analysis  has  been  applied  to  a variety  oi  infrared  mirror  and  window  surfaces  pre- 
pared by  various  techniques.  Both  coated  and  'mcoated  samples  have  been  studied.  The  present  section 
gives  a few  representative  examples  of  results  obtained. 

The  application  of  multithreshold  analysis  tc  relate  damage  mechanisms  to  surface  characteristics 
of  bare  metal  mirrors,  particularly  Cu,  has  beer,  amply  demonstrated  [lj.  Certain  aspects  of  the  Au 
mirror  represented  in  figure  3 were  also  discussed.-  To  complete  this  discussion  it  should  be  pointed 
out  that  the  relatively  high  slip  threshold  indicated  in  figure  3 may  be  influenced  by  the  difficulty  in 
observing  slip  on  this  sample  [1],  In  fact,  the  lowest  threshold  is  for  work-function  change,  which  is 
attributed  to  evolution  of  water  by  the  Au  layer..  This  result  is  supported  by  a laser-induced  diminu- 
tion of  oxygen  as  observed  by  Auger. 

The  value  of  multithreshold  analysis  in  evaluating  nominally  eouivalent  mirror  finishes  is  illus- 
trated in  figure  4,  where  we  compare  damage  profiles  of  Mo  mirrors  polished  by  two  different  vendors. 

While  mirror  A exhibits  higher  thresholds  for  most  types  of  damsge,  including  melting,  it  is  far  more 
susceptible  to  pitting  and  associated  plasma  formation.  The  superior  uniformity  of  mirror  B is  also 
evident  in  the  smaller  standard  deviations. 

As  an  example  of  multithreshold  results  from  a window  material  we  compare  profiles  for  damage 
occurring  on  the  entrance  (front)  and  exit  (back)  surfaces  of  3 polished  NaCl  blank  (fig.  5).  Only 
chose  sites  where  no  simultaneous  damage  was  observed  on  the  front  surface  are  represented  in  figure  5(b).; 
All  thresholds  on  the  front  are  essentially  the  same.  Indicating  breakdown  as  the  only  damage 
mechanism.  Thresholds  at  the  back  surface  are  generally  lower  because  of  the  constructive  phase  rela- 
tionship between  incident  and  reflected  fields.  Effects  peculiar  to  Che  back  surface  are  ripples,  pits 
snd  fracture.  Ripples  result  from  the  interference  of  incident  and  scattered  radiation  [6].  Present 
evidence  shows  that  ripples  occur  only  in  the  presence  of  a plasma.  The  pits  may  represent  scattering 
sources  for  the  ripples,  '"he  fact  chat  pita  occur  only  on  the  back  surface  suggest  that  they  result  " 
from  bulk,  ratner  than  jb  surface  inbonogeneities,  Scrubbing  associated  with  the  plasma,  vhiuh  was 
observed  in  mirrors  [1],  has  also  been  observed  on  this  saaole.; 

A damage  profile  from  an  enhancad-reflectiun  coating  on  polished  bulk  Cu  is  shown  in  figure  6.-  The 
thresholds  are  well  below  those  for  bare  Cu,  Pitrirg  cue  to  coatirn,  defects  is  the  dominant  failure  mode-. 
Almost  identical  damage  piorilcs  were  obtained  when  the  same  coatm;  design  was  applied  to  Mo  and  to  ULE 
glass  .- 

Figure  7 compares  damage  thresholds  ffca  as^Sj  in  buck  ferr.,  and  as  a costing  (?]  on  the  NaC'l.  In 
the  bulk  material  the  coincidence  of  thrashoids  for  surface  pitting  and  bulk  damage  suggest  that  these 
result  from  the  same  mechanism  which,  surprisingly,  is  active  below  the  threshold  for  pitting  m the 
coating.  Also,  pitting  in  the  coating  tends  to  associate  with  a plasma,  wnich  in  nor  Lrue  in  the  bulk 
material.; 


4 .  Summary 

An  efficient  and  reliable  method  has  been  developed  for  simultaneously  measuring  the  thresh- 
olds and  standard  deviations  for  nr>  to  eight  laaer-damage-reieted  srfecto  on  a single  sample.  An 
algorithm  based  on  the  maximum-likelihood  principle  has  beer,  derived  for  computing  results  directly  from 
the  experimental  data.  Reproduclb '.ilry  has  eoer.  demonstrated  by  repeated  measurements  on  a standard  Au 
sample.  The  use  of  this  sample  to  maintain  long-term  calibration  was  described.-  Txamples  of  multi- 
threshold results  on  representative  optical  component-  of  each  tvoe  haw  been  given  and  the  principal 
weaknesses  of  these  components  identified.  Tbe  comparative  analysis  o*  such  results  Cron  well  < harac- 
terized  samples  prepared  from  different  materials  and  bv  diffeient  techniques  yields  the  relationship 
between  damage  thresholds  and  material  cherectei  istics , These  compprloep.s  can  also  provide  a better 
understanding  of  laser-damage  mechanisms. 
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7.  Figures 
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Figure  1 . Micrographs 
of  visible  damage  on 
polished  bulk  Cu  at  a 
peak  fluence  ol  68.3 
J/cm2.  In  Nomarski, 
slip  bands  appear 
outside  the  crater  at 
the  5 o'clock  position, 
while  intergranular 
slip  is  responsible  for 
prominent  crack-like 
features . 


Figure  2.  Gaussian  spatial 
profile  of  focussed  laser 
pulse,  showing  1/e2  width  and 
peak  energy  density  (fluence) 
at  which  thresholds  are 
determined . 
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Figure  3(a).  Damage  profile  of  standard  gold  sample, 
measured  for  calibration  with  the  original  laser  out- 
put coupler.  Thresholds  are  indicated  by  the  shaded 
bars.  Standard  deviations  (■*■  or  -)  are  indicated  by 
the  semishaded  boxes  at  the  ends  of  the  bars. 
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Figure  3(b).  Damage  profile  to  be  compared  with 
figure  3(a),  measured  for  calibration  with  new 
laser  output  coupler. 
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Figure  6.  Damage  profile  of  a 10.6-ym  enhanced  reflection 
coatir?  on  bulk  Cu.  Coating  consists  of  alternating  layers 
of  Ge  and  ZnS  on  a thin  Ag  layer. 
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Figure  7(a).  Damage  profile  of  an  AS2S3  coating  on 
the  front  surface  of  NaCl. 
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Figure  7(b) . 
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